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ABSTRACT 

We introduce a new technique to quantify highly structured spectra for which 
the definition of continua or spectral features in the observed flux spectra is dif- 
ficult. The method employs wavelet transformation which allows the decompo- 
sition of the observed spectra into different scales. A procedure is formulated to 
define the strength of spectral features so that the measured spectral indices are 
independent of the flux levels and are insensitive to the definition of continuum 
and also to reddening. This technique is applied to Type la supernovae spectra, 
where correlations are revealed between the luminosity and spectral features. The 
current technique may allow for luminosity corrections based on spectral features 
in the use of Type la supernovae as cosmological probe. 

Subject headings: Cosmology: Distance Scales, Stars: Supernova, Methods: Ob- 
servational, Techniques: Spectroscopic 



Introduction 



A major difficulty in analyzing spectroscopic data with highly blended atomic lines is to 
quantify the strength of certain spectral features. These spectral features are superimposed 
on a continuum so line blending can make it difficult to reliably define the continuum level. 
In the case of supernova spectra, the pre-nebular phase spectra typically show P-Cygni 
profiles with both emission and absorption components whereas the nebular phase spectra 
are dominated by broad overlapping emission lines. The spectral features are therefore of 
various widths and strengths, and neighboring features are heavily blended. Further, the 
data typically contain observational noise, flux calibrations errors and uncertainties in the 
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amount of dust extinction. The noise makes the definition of a continuum very uncertain and 
accordingly the calculation of equivalent width becomes unreliable. In many observations, in 
particular those at high redshift, the observed supernova spectra are heavily contaminated 
by host galaxy spectra. This affects severely the definition of line depth. 

For Type la supernovae, it is known that certain spectral line ratios such as the Si 
II 5800/Si II 6150, and the Ca I I H&K lines are sensitive to the intrinsic brightness of 



the supernova iNugent et al.l (119951 ). The measurement of the line strength is, however, not 



trivial . For instance, to measure the Si II 5800/Si II 6150, and the Ca II ratio, INugent et al. 



(119951 ) employed a simple approach by drawing straight lines at the local peaks of the spectral 
features and measure the depth of the absorption minima from the straight line. However, 
the location of the straight line and the position of the line minimum are not easy to define 
in the presence of observational errors. It is for this reason the ratio is derived only for a 
number of very well observed local supernovae. 

In this paper, the spectral features of Type la supernovae will be analyzed instead 
through wavelet transformations. This technique avoids many of the challenges mentioned 
above associated with identifying line strengths. Here wavelet transformations are applied 
to Type la supernovae spectra with the purpose to quantify the spectral features for cosmo- 
logical applications. 



2. Spectral Features 
2.1. Wavelet Transform Algorithm 



It has been pre v iously demo nstrated that the a trous algorithm iHoltschneider et al. 



(1989 


); 


Starck et al. 


(1995. 


1997) 



transform is carried out in direct space so artifacts related to periodicity do not occur. The 
reconstruction is trivial. The evolution of the transform from one scale to the next is easy 
to follow and the interpretation of the spectrum at each scale is straightforward. 

The a trous wavelet uses a dyadic wavelet to merge non-dyadic data in a simple and 
efficient procedure. Assuming a scaling function <p(x) (which corresponds to a low pass filter), 
the first filtering is performed on the original data {co(k)} by a twice magnified scale leading 
to the {ci(k)} set. The signal difference {co(k)} — {c\{k)} contains the information between 
these two scales and is the discrete set associated with the wavelet transform corresponding 
to 4>(x). The operation is performed successively and to obtain the wavelet scale Wj{k) at 
each scale j. The original spectrum cq can be expressed as the sum of all the wavelet scales 
and the last smoothed array c p : 
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Fig. 1. — (a) The spectrum of SN 2001el, (b-j) the wavelet scales 1 to 9, and (k) the smoothed 
array c p of the SN 2001el from top left to lower right. A sum of of scales (b) to (k) recovers the 
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co(k) = Cp(k) + TP j=1 Wj(k) 

To demonstrate the basic features of the a trous wavelet transformation, we show in 
Figure [1] the wavelet transformation of a well observed supernova SN 2001el. The data 
were obtained through the spec tropolarimet r y pro gram at the Very Large Telescope of the 



European Southern Observatory I Wang et al.l (120031 ). The sampling step of the data is binned 



to 5A. The signal to noise ratio (SNR) of the data is everywhere above 150 - this unusually 
high SNR is a result of the spectropolarimetry observations. The original data is show in 
Figure [Jji, and the consecutive wavelet scales for j = 1, 9 are shown in Figure [lb to [Tj. 
Figure [Ik represents the smoothed array c p : given the spectral range of SN 200 lei the a 
trous wavelet cannot generate more than 10 wavelet scales. Further note that each of the 
individual wavelet scales have zero mean. It can be seen that at small scales the wavelet 
is dominated by observational noise and the supernova signal starts to become significant 
only for j > 3, and the broad spectral wiggles associated with the supernova dominate the 
wavelet scales of j = 5, 6, and 7. The supernova spectral features are typically a few hundred 
A wide and are effectively isolated in the decomposed spectra. 

The spectral features of a supernova can be better described by a blend of several 
wavelet scales. For this reason, we can calculate the sum of more than one scales to reflect 
the existence of features of various width: 

W {1} = ZjeyyWj, 

where {/} is a subset of wavelet scales. Examples of these spectra are shown in Figure [2] for 
SN 2001el. 



2.2. Normalization of Spectral Features 

The wavelet scales, having the units of the original flux spectrum, need to be normalized 
to construct quantities that measure the strength of the spectral features that do not depend 
on the absolute flux level of the spectrum. There undoubtedly is more than one way to 
normalize the scales. The simplest approach is to normalize all the wavelet scales by dividing 
them by the smoothed array c p . This approach is simple and will certainly work fine for data 
without host galaxy contamination. For data with host galaxy contamination, or those with 
poor background subtraction, this approach introduces systematic errors to the normalized 
scales. 

In our approach, the normalized wavelet scale is defined using the standard deviations 
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Fig. 2.— Wavelet sums for SN 2001el over (a) scales 1 to 3, (b) 2 to 4, (c) 3 to 5, (d) 4 to 
6, (e) 5 to 7, (f) 6 to 8 and (g) 7 to 9. Running combinations of the different scales captures 
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of the spectral features from any given wavelet scale: 

W m (X) = W {l} (\)/^Wf l} (\)/N l2 (1) 
= W {l} (X)/a {l} , (2) 

where N i2 is the number of data points between Ai and A 2 . The mean and standard deviation 
of W{iy are zero and 1, respectively. This is effectively a self- normalization that exploits only 
the intrinsic properties of the wavelet scales involved. Host galaxy contamination (which 
does strongly affect c p ) would not have a significant effect in this context. 

The spectral index Xj of any feature between A a and \ b at a given scale j is defined by 
averaging the normalized wavelet scale Wf 

X {1} = ^»W {l} (\)/N ab , 

with N ab being the bin size in the wavelength region A a and X b . Xj defines a normalized 
number which measures the strength of the spectral features in the normalized wavelet scale 
Wj between A a and X b . Alternatively, one can also calculate the power Pj between A a and 
X b for wavelet scale j: 

p {i} = ^W 2 {l} (X)/N ab . 
P{i} and X{;} contain the same information. In this study, we will focus on Xyy. 

One obvious advantage to using wavelet scales to estimate spectral feature strengths is 
that they do not depend on the definition of the spectral continua. Furthermore, since they 
can be estimated locally around a spectral feature, spectral indices are useful in minimizing 
uncertainties due to errors in spectral flux calibrations. Similarly, the spectral indices as 
defined here are less sensitive to errors of background subtraction, which is usually one of 
the dominate sources of uncertainty, especially in the studies of high redshift supernovae. 



2.3. Normalization Spectral Features of SN la 

The wavelet technique is particularly well-suited for studying scattering-dominated spec- 
tra of expanding atmospheres with P-Cygni spectral features: the net flux of the P-Cygni 
feature is usually close to zero. Wavelet decomposition is consistent with this as the mean 
flux is zero for the various wavelet scales. Wavelet transforms thus makes it easy to separate 
emission and absorption components of a spectrum in a mathematically robust way. 

In this study, the supernova spectra are first decomposed into various scales as described 
in the above section. In addition, to reflect the fact the spectral features are a blend of 
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different scales, the sum of the wavelet scales 3, 4, and 5 are used as the primary spectrum 
for the analysis of spectral features (though other scales have also been analyzed). All 
the decomposed spectra are normalized in a similar way as given in Equation [2J To derive 
quantities that are less sensitive to errors of flux calibration, we need to restrict calculation of 
the normalization factor to a small wavelength region and yet to have large enough spectral 
coverage so that the feature strengths will not be affected by boundary. In this study, the 
spectra are divided into four regions: (A) 5500 to 6500 A, (B) 4985 to 5985 A, (C) 4850 to 
5450 A, and (D) 4250 to 5200 A; the variance in each of these sections of spectra is calculated 
and used as the normalization factor. Interesting features include the Silicon II lines at 635.5 
and 580.0 nm in region (A), the Silicon II line around 548. 5nm in region (B), and the two 
strong peaks at 510.0 nm and 450.0 nm, in region (C) and (D), respectively. These five 
spectral features are shown in Figure [3] and are the main focus of subsequent analyses. 

3. Biases and Errors 

In practice, the observed data contain noise and estimates of Xj can be biased. The 
noise affects Xj in two ways: First, when the noise is large, its effect can propagate to 
all the wavelet scales and become a significant component at the wavelet scale of interest. 
Secondly, it changes the normalization factor when calculating Wj - data with larger noise 
can be systematically biased to give a larger normalization factor because the additional 
power from shot noise. This bias is usually not a problem for high SNR ratio data, but can 
be significant for data with a low SNR ratio. The correction factor Tl(j) for scale j is defined 
as: 

a (j) =a(j)U lv (3) 

where cr (j) is the variance at the jth scale in the ideal case of no photon shot noise. 

Typically, published spectroscopic data do not have associated noise spectra. One in- 
stead has to rely on the flux spectrum to estimate the noise levels. A major advantage 
of wavelet transformation is that it allows estimates of the noise characteristics based on 
the spectral data itself. If we assume that all the continuum or spectral features are much 
broader than the data sampling step, the spectral fluctuations of the wavelet scale with 
j = 1 should then represent mostly shot noise. This is generally reasonable as can be seen 
in Figured] (b) for the spectrum of Type la supernova 2001el. 

Recognizing that smaller wavelet scales contain more information of the noise property 
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Fig. 3. — Spectral features for Type la supernovae. (a): 1 = 1,2,3. (b): 1 = 2,3,4, (c) 1 
3,4,5, and (d) 1 = 4,5,6. The spectral features are well resolved in (c) and (d). 
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than larger wavelet scales, we can define the spectral quality index (SQI) of the normalized 
wavelet scale {/} as the variance ratio of normalized scale {1} and the lowest normalized 
scale {1}: 

p {l}{1} = T^Wy^W^, (4) 

where the braces {} reflect that the various quantities are actually sums of various 
wavelet scales. Specifically, W{i} is the normalized sum of three wavelet scales, where {1} 
= 1,2,3 in this instance. The SQI measures the relative importance of noise levels in 
estimating of the spectral feature index. It can be calculated directly from the decomposed 
spectra without an error spectrum. Note that SQI is a quantity that can be localized to 
certain wavelength intervals. 

For a given spectrum {q}, the dependence of the wavelet spectral indices Xyy and the 
correction factor II on SQI can be estimated through Monte-Carlo simulations. 

3.1. Dependence of Spectral Features on Observational Noise 

Monte Carlo simulations are required to quantify the dependence of X indices on obser- 
vational noise. The characteristic parameter of the noise is the SQI defined in Equation [3]- 
the ratio of spectral variance of combined wavelet scales I = 3, 4, 5 to that of the combined 
scales I = 1,2, 3. To perform such these simulations one needs a series of noise free spec- 
tra of supernova spectra. The spectropolarimetry program at the ESO VLT has acquired 
several high quality spectra of SN la with SNR ratios around 150 (2003ApJ...591.1110W). 
Spectra of SN 2001V and SN 2001el from the spectropolarimetry program will be used in 
this simulation to quantify the relations between X and SQI. 

In the example shown in Figure HI various levels of Poisson noises were added to the 
spectrum of SN 2001el at day +1. The noise is added to the spectra, which are then 
transformed to various wavelet scales and the various X indices are calculated. The top 
panel in Figure H] shows the relation between p and the assumed SNR ratio with the addition 
of Poisson noise. The SQI is calculated in the wavelength intervals of 550.0 nm to 650.0 nm, 
498. 5nm to 598.5 nm, and 425.0 to 520.0 nm. It can be seen that p is correlates well with 
the SNR ratio of the input data: reducing the SNR ratio decreases p. This confirms that 
the SQI can effectively capture effect the photon shot noise, and can be used to quantify the 
noise level of the data. 
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The variances used to normalized the spectra at the various wavelet scales are clearly 
correlated. This is shown in the middle panel of Figure HI where the data exhibit nearly 
identical slopes over the different wavelength regions. A linear relationship between cx 2 (l) 
and cr 2 (3), and between a 2 (l) and cr 2 (4) is assumed for the fits. The slopes 7^ extracted 
from these fits are given in Table [TJ 

The bottom panel in Figure H] clearly demonstrates how the correction factor II increases 
dramatically for p approaching 2.82 (which corresponds to a SNR ratio of below 1 per 0.5 nm 
bin). This implies that the spectral features are dominated by the noise, hence it becomes 
impossible to extract the spectral indices reliably. 

The correction factor for p can be fit well with a function 

n v = (5) 

with the relevant coefficients taken from Table [2] for the various lines. 



3.2. Bias corrections 

The various X indices for the spectral features are derived from the Monte-Carlo sim- 
ulation of data with different SQL As shown in Figure El the X indices (shown as open 
squares apparently suffer strong bias when the data are noisy. The various X indices after IT 
corrections are shown in Figure [51 The effect is generally small for high SNR ratio data, but 
becomes important for data with low SNR ratio. In any case the bias is effectively removed 
by applying the correction factor IT. 



3.3. Error Estimates of the Spectral Indices 

Assuming photon shot noise, the Monte Carlo simulations also give error estimates for 
the X indices. The errors as a function of p are shown in Figure [6j These errors are fitted 
with a function of the form: 

ox = (6) 



and the relevant coefficients 77 and ip are shown in Table [3 Simulations were performed 
for all of the SN 2001V and SN 200 lei spectra and it was found that in all cases the bias 
can be well corrected. Note that due to the lack of a completely noise-free SN la spectrum, 
at extremely high SNR ratio (such as those that are higher than or comparable with the 



- 11 - 



5500-6500 



-1985-5985 



4850-5450 



4250-5200 



10 100 1 10 100 1 10 100 1 10 100 

S/N Ratio S/N Ratio S/N Rotio S/N Rotio 





/ 


/ - 


y 

/ 



10 20 30 40 50 



10 20 30 40 50 SB 10 20 30 40 50 SB 10 20 30 40 50 60 
a 2 (l) a 2 (l) a 2 (l) 



12 3 4 5 
1 r m 





234501 23450 
VP, 1/P M 



2 3 4 5 
'//>„ 



Fig. 4. — (Top) The relation between the SQI and the input SNR ratios. The solid lines show 
SQI of fourth wavelet scale and the dashed line the third. (Middle) The relation between the 
variance of the third and first wavelet scales. The effect of a large er(l) propagates linearly 
to larger wavelet scales. (Bottom) The correlation of the bias correction factor and SQI for 
the third wavelet scale. The SNR was varied in all cases via the addition of noise in the 
Monte Carlo simulations. 
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Fig. 5. — Line indices corrected for p dependence for important SN Type la spectral features. 
The X indices are derived from the sum of of wavelets 3,4,5 of SN 2001el 1 day past optical 
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signal to noise ratio of the SN 2001el spectra as used in the simulation) the Monte Carlo 
simulations do not give correct estimates of the errors. Such cases are unlikely to be relevant 
as in such situations, the errors are likely to be dominated by calibration systematics rather 
than shot noise. The error function given above will be used for all cases here. As can be 
seen in Figure El the above expression gives an excellent description of the dependence when 
the errors are described by p. 

The p dependence of the X indices and their errors have a weak dependence among the 
different varieties of supernovae and the epoch of the supernovae. 

Procedure for Bias and Error Estimation 

The procedure for removing bias and estimating errors from noisy supernovae spectra is 
premised on extracting correction factors from a supernova with a large SNR. Here we 
enumerate a correction recipe using wavelet scales I = 3, 4, 5 with SN 2001 el as our reference 
spectrum. 



1. Compute the ratio of the sum of the squares of wavelet scales 3,4,5 and 1,2,3 (p 3i in 
Eq. E} for SN 2001 el 

Km 

P31 = -r 1 -• 

2. Degrade the SNR ratio of SN 2001 el in multiple steps with the addition of Poisson 
noise. Compute aj x 2 3 i and <7| 3 4 5 | at each step. 

3. Extract 731 by assuming a linear relationship between cT{ 12 3} an d o"{ 34 5 } (see middle 
panel of Figure 0J. 



^{3,4,5} = /5 + 731 * cr {l,2,3}> 

where a^} is defined in Eq. |2j 

4. Repeat each of the above steps over all regions of interest to extract a mean value for 
731. Equivalently, use the values for 7 31 by consulting Table [TJ 



Table 1: The Coefficients for the Dependence of X on Data Errors 



SN Day 


731(A) 731 (B) 731(C) 731(D) 


mean 


741(A) 741 (B) 741(C) 741(D) 


mean 


01V -8 
Olel -4 
Olel +1 
Olel +9 


0.1132 0.1116 0.1183 0.1144 
0.1220 0.1236 0.1229 0.1224 
0.1220 0.1236 0.1229 0.1224 
0.1251 0.1273 0.1250 0.1240 


0.1144 
0.1227 
0.1227 
0.1254 


0.0263 0.0240 0.0273 0.0259 
0.0567 0.0587 0.0566 0.0573 
0.0581 0.0583 0.0552 0.0551 
0.0585 0.0607 0.0599 0.0577 


0.0258 
0.0573 
0.0567 
0.0592 
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Fig. 6. — Errors of the X indices as functions of p for SN 2001el at optical maximum. 
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5. Using equation 13. 11 compute p$i values for each SNe having typical values of the SNR 
ratio. 

6. Compute the correction factor n 31 using this value of p 3i and the value of 731 computed 
for SN 2001 el 



n 3 i = yl-TsiPsi- 

7. The spectral index of any supernovae feature can be corrected for bias by dividing the 
uncorrected value by II31 as in Eq. 

x x 



tt31 

8. The error bars are determined from the same set of simulations. With ox defined as 
in Equation 13.31 construct the relationship 

log w (a x ) =V + iplog w (p 31 ), 

where all quantities refer to a SNe with a large SNR ratio (e.g. SN 2001 el). Equiva- 
lently, use the values for 7 31 by consulting Table [2j 

9. With these values i] and ip compute the variance in the spectral index for a SN with a 
typical value of SNR as 

a x = lOVt 



4. Applications to Type la Supernovae 



Spectral indices lend themselves to a quantita t ive an alysis of the temporal and magni- 
tude evolution of the spectral lines. iNugent et al.l (119951 ) measured the ratio of the depths 
of the Si II 6150 A and 5750 A features and established correlations with Amis. Other 
studies of these and other spectral fea t ures h ave adopted slightly more elaborate p rocedures 



based on equivalent lHachinger et al.l (12006!) and p s eudo equivalent wid ths (EW) iFolatelli 



(|2003h : lGaravini et al.l (ho07i ): lAltavirla et all (l2009h: 



Branch et al 



J2009h to stu d y the se and 



other absorption features. Recently, lArsenijevic et al. ( 20081 ) ; Stanishev et al. ( 2007 ) have 



Table 2: The Coefficients for the Errors o 



X 



SN 


Date 


??(6150) 


i/>(6150) 


>7(5800) 


i/>(5800) 


?7(5485) 


■0(5485) 


>7(5150) 


■0(5150) 


>7(4250) 


i/>(4250) 


2001V 


-8 


-0.02486 


0.0469 


-0.0156 


0.0293 


-0.0203 


0.0601 


-0.0373 


0.0545 


-0.0179 


0.0453 


2001el 


-4 


-0.00625 


0. 0367 


-0.00628 


0.05429 


-0.02449 


0.05892 


-0.01626 


0.06382 


-0.01165 


0.05248 


2001el 


+ 1 


0.133 


1.248 


0.212 


1.188 


0.160 


1.412 


0.176 


1.533 


0.157 


1.531 


2001cl 


+9 


0.159 


1.276 


0.194 


1.282 


0.234 


1.279 


0.203 


1.442 


0.177 


1.672 


2001cl 


+ 1 


0.145 


1.194 


0.265 


1.322 


0.175 


1.295 


0.181 


1 . 500 


0.216 


1.408 
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used wavelets coupled with the pseudo-EW technique to study a Si II absorption feature. 
Emission features have received less attention t han absorption featur e s and usually involve a 
distinct p rocedure from the a bsorption features iNugent et al.l (119951 ); iBongard et al.l (120061 ). 
Recently, iBailey et al.l (120091 ) described a variance of the above methods wherein absorption 
and emission features in a training set of spectra are studied to extract the optimal flux 
ratio to Amis correlation. This ratio is then applied to correct the magnitudes of other 
supernovae within a validation set. 

The wavelet technique developed here differs in several important respects from those 
described in the preceding paragraph. Our methodology is premised on the existence of one 
(or more) very high SNR spectra. Spectral line strengths are first extracted for this high 
SNR spectrum from a combination of intermediate wavelet scales. Excluding the lowest and 
highest reduces the effects of noise and the continuum, respectively. Working with wavelets 
from a high SNR spectrum allows corrections to be made to lower SNR spectra. Perhaps the 
most salient difference is that performing our analysis entirely in wavelet space permits us 
to avoid definition of the continuum (the mean of the wavelet scales is zero and integration 
is performed from one zero on the leading to the next zero on the trailing edge of a feature). 
Consequently, we are able to work directly with line strengths of the features themselves, 
not their ratios. Lastly, this technique permits absorption and emission lines to be treated 
democratically. It is our expectation that the wavelet method gives to a more robust measure 
of line strength. In Sections 14.21 and 14.31 we apply this technique to data described in the 
subsequent section. 



4.1. Data Sample 

The supernovae included in this study are given in Table 14.11 A large number of low- 
z spectra (z < 0.1) were collected from libraries that are publicly available, such as the 
SUSPECT Supernova Databas^j] and the Center for Astrophysics Supernova Archive^], as well 
as other SNe that are available in the literature. The spectra are corrected by the host galaxy 
redshift but no dust extinction correction is applied. The original wavelength coverages, step 
sizes and SNR ratios of these spectra are vastly different. In our analysis, all the spectra 
are first rebinned to 5 A sampling step for convenience. After wavelet decomposition was 
performed the spectra were checked for edge effects that would distort calculations, the 
affected spectra were removed. 



1 http://bruford. nhn.ou.edu/ suspect/indexl.html 
2 http://www. cfa.harvard.edu/supernova/SNarchive. html 
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Table 3. The Spectroscopic Sample of SNe la 



SNc 


Ami 5 a 


Branch Subtype* 3 


^6150 


^5750 


-^5485 


^5150 


X4570 


h])(< 1 1 a Source 


1981B 


1.125(0.010) 


BL 


-56.645(0.157) 


-10.339(0.298) 


-48.941(0.245) 


42.956(0.175) 


43.306(0.184) 


1 


1983 G 


1.37{0.01j b 




-59.818(0.649) 


-13.890(3.411) 


k 


49.532(0.483) 


k 


2 


1984A 


1.294(0.063) 


BL 


-63.615(0.974) 


-5.166(1.167) 


-39.319(0.489) 


45.827(0.573) 


43.895(0.781) 


3 


1986G 


1.643(0.022) 


CL 


-49.109(0.295) 


-20.165(0.534) 


-32.227(0.456) 


40.231(0.386) 


51.494(0.521) 


4 


1989B 


1.262(0.017) 


CL 


-54.202(0.486) 


-13.218(0.428) 


-47.078(0.424) 
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a Ami5 values were calculated by the super-stretch method from [Wang et al ] 120061 ) unless otherwise noted, 
k Designations from [Branch et al.| j2009|) 



- 18 - 



c Ami5 from [Hachinger et aT ] 120061 1 
d Ami 5 from lJha et al.l 11999ft 
°Am 15 from iPhillips et all 120071) 
' Amu from lPignata et all 120081) 

g Ami5 converted from stretch value, s, from Quimby et al. (2006) using the equation from Pcrlmuttcr et al 
Amis converted from stretch value, s, from Quimby ct al. (2007) using the equation from Pcrlmuttcr ct al 
'Amis from Hickcn ct al ] 120071) 
j Am 15 from I Wang et al. (2008) 

k The 5485 A and 4570 A features show much more variance in their evolution, therefore the epoch range over which these features were fit was smaller. 
These SNc arc missing ^5485 and -^4570 values because they did not have enough spectra within the smaller epoch range. 

'Due to noise or miscalibration of the spectra at +3 days, there is not enough data to fit -Y5485, -^5150' and ^4570 

m The spectra for SN 2001el did not cover the wavelength region for this feature. 

n Not enough of the spectra for SN 2005bl covered the wavelength regions for X5485 , -^5150, and ^4570 for a good fit to be made. 



Ref eren ces. — fl) | Branch et aP 119831); (2HHarriB et all ft983l ); (3)lBarbon et aljll989l) ; (4) lPhillips et al.l tl987l); )5) [Barbon et all 119901); (6)IWells et~al"l 
119941); f7)IMazzali et al.l 11993T); f8)]Leiburidgut et al.1 119911) ; f9)lFilippenko et al.[fl992a[); (lO)lPhillips et al.lll992i) ; (ll) lRuiz-Lapuente et allll992|); (12) 
iLeibundgut et al.l 119931) ; (13) |Firippenko et alj 11992bl) ; "( 14) iTuratto et al.l 119961); f 15) ItCh-sImer et all 119931) ; (16 ) Wang unpubli shed; ( 17) IMeikle e~ 
199c k (18) iPatat et al.l tl996h ; ( 19) [Wang et al l 119971) ; ( 20) ISalvo et al I 12001ft: (21 ) lAnupamal 11997ft; (22) iLi et al.l ll999l); (23) iMatheson et . 
200S ): (24) IBranch et al.l 12003): (25) Ijha et al I ll999lj: (26) IMeikle fc Hernandez! 120001): (27) iHernandez et alj 120001) ; (28) 1 IGaravini et al.l l2004fl; 
IGaravini et alj 120051); f 30) IStrolger et alj 12003) ; (31) iHamuv et al.l 120021) ; (32) IValentini et al.1 I2003D ; (33) ILi et alT 120011) ; (34) IBranch et al.|"lf2"006l 
(35)IWang et al.1 120031) ; (36 ) IBenetti et al J 120041) ; (37) ILi et al.l 120031); (38) | Pignata et alj 12008); (39) Wan g unpublished; (40) iKotak et al. | |2005l); (41) 
Elias-Rosa ct al. (2006); (42) Anupama ct al. (2005); (43) Stanishcv ct al. 12007); (44) Howell ct al. 12006); (45) Krisciunas ct al. 12007); (46) Altavilla ct al. 
120071); (471 iMazzali et al-l 120081); (481 iTaubenberger et al.l 12008) ) ; (49) IGaravini et alj 12007|) ; (50) IQuimbv et al.1 120061) ; (51) Quain in progress; (52) 
IQuimbv et alj 120071) ; (53) Iphillips et al.1 120071 ); (54) lHicken et all 120071 ); (55) IWang et al.1 
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4.2. X Versus the Epochs 

An example of the time evolution of X indices is shown in Figure [7] for SN 2005cf - a 
normal Type la supernova with Amis — 1.16. For a spectroscopically normal supernovae like 
SN 2005cf, the Si II 5800 and 6150 lines exhibit little evolution in line strength for roughly 
±8 days around maximum, after which the the 5800 A line strengthens and the 6150 A line 
weakens. Similarly, around 8 days past maximum the S II feature begins to weaken until it 
is completely obscured by 18 days past maximum. The emission features at 4750 and 5150 
A for this same supernova, by contrast, shows comparatively little time evolution. 

Analysis of the time evolution is complicated by occasional large gaps between epochs 
and the need to occasionally track spectral features manually due to the decreasing velocity 
of the expanding photosphere. Consequently, a full analysis of the temporal evolution of the 
remaining supernovae in Table 14.11 will be analyzed in a separate paper. 



4.3. X versus Ami 5 

Figures 181 to [T21 show the correlations between X and Am 15 for the five spectral features 
that we have adopted for our study. The line strengths in these figures are those computed at 
maximum light. In instances where no spectrum at maximum light exists a simple quadratic 
fit was made of all spectra within 8 days of maximum (for features that are not a smoothly 
varying the fit was restricted to within 5 days of maximum). The fits were checked for 
consistency and for supernovae with only two spectra closely sampled in time, a mean was 
taken to avoid aberrant behavior in the fit. Any supernova having only a single spectrum 
within the specified time range was removed, unless that spectrum was taken at maximum. 



4.3.1. The Si II 6150 A line 



It has been shown previously that the strength of t he Si II 6150 A line is not tightly 
correlated with the intrinsic brightness lHachinger et all (120061 ) . Figure [8] confirms this ob- 
servation in the main. However, the X indices do show a modest trend of weaker spectral 
strength for dimmer supernovae and supernovae with Amis less than 1 show large variations 
of the Si II 6150 strength. It merits mention that the several of these supernovae (e.g., 
SN 1997br, SN 200 lay, SN 2002cx, and SN 2005hk) are deviant with respect to the majority 
of the sample. It has been noted that these are all peculiar supernovae and it has been 
speculated that SN 1997br, SN 2002cx , and S N 200 5 hk may form a group distinct from most 
typical Type la supernovae iLi et all (ll999l . 120031 ); iBranch et all ( 120041 ); iHowell fc Nugent 
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Evolution of Spectral Indexes for SN 2005cf 
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Fig. 7. — The temporal evolution of the X indices of SN 2005cf. In order these are: (a) the 
emission feature at 4570, (b) the emission feature at 5150, (c) S II line, (d)Si II 5800 line, 
(e) the Si II 6150. 
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Fig. 8. — The correlation of the strength of Si II 6150 feature and Am^. 
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Branch et al 



(I2004h: Branch et all (120061 ): Ijha et al.l (J200J); IPhillips et all (120071 ): ISahu et all (120081 ): 
(l2009h . 



^.3.2. Si II 5800 A 

It has previously been shown that the ratio of t he strength of this fe a ture and that of 



the Si II 6150 A line are well correlated with Amis iNugent et al.l (119951 ): iHachinger et al. 
(120061 ). Using this correlation, a determination of t he supernoy a es' m aximum luminosity 
may be determined on the basis of a single spectra iRiess et al.l (119981) . These same two 



features hav e also been used to define la SNe subgroups iBenetti et al.l (120051 ): iBranch et al. 

(boodboogh . 



Figure [9] shows the correlations with Amis. The X index defined here for this line 
correlates tightly with Am 15 even without having been divided by the strength of the Si II 
6150 line. Note however, that the X index for this feature is normalized by the total variance 
of the wavelength scales from wavelength region between 5500 A and 6500 A, the variations 
due to Si II 6150 A line are partially included in the derivations of the X indices. The 
correlation between X and Amis can be well described by a linear relation. It may be the 
case that for sub-luminous s upernovae such a s SN 1 991bg, the feature 5800 Amay actually 
be a blend of Ti II a nd Si IllGarnavich et al.l (12004 ) (although this has also been disputed 
Branch et all (|2006h : feongard et all (bo08h ). Blending would artificially enhance the line 
strength and may suggest more of a correlation than in the absence of line blending. We 
wish to emphasize that the X index of this feature measures the total strength of this feature 
and does not distinguish the physical origins of the features. 



4.3.3. S II 5485 

The "w" shaped spectral f eature S II 5485 i s ano ther important line that defines an SN la 
similar to the Si II 6150 lines iBongard et al.l (120061 ) . Figure [TUl suggests that the strength 
of this feature too may depend on Am 15 , though correlation is much weaker than that for 
Si II 5800. The deviant SNe are SN 1991T, SN 2001ay, SN 2005hk, SN 2006X, SN 2001V 
and SN 2002cx. SN 1997br has a similar spectral evolution to that of SN 1991T, but it has 
not been shown since it has only one spectrum within 5 days of maximum However the X 
value (—30) for SN 1997br at 4 days before maximum is consistent with the X value for 
SN 1991T. Other 1991T-like SNe (SN 1998es, SN 1999aa, SN 1999dq, SN 2000cx) have more 
normal values but they still are on the upper edge of the distribution. There is apparently 
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Fig. 9. — The correlation of the strength of Si II 5800 and A77115. 
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Fig. 10. — The correlation of the strength of S II 5485 and Amis. 
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some diversity within the so-called SN 1991T-like SNe. 

The S II line is generally much stronger than the Si II 5800 line and is thus much easier 
to measure. There is also much more evolution within this feature so Figure 10 is restricted 
to spectra taken within 5 days of maximum. 

4.4. 5150 Emission Feature 

Figure [TT1 shows what might be interpreted as a slight dependence on Amis, though the 
trend is not as apparent as it is with the S II line. The X values do decrease notably with 
Amis for the most sub-luminous supernovae. A majority of the data are clustered with a 
large amount of scatter. The 5150 A emission feature may not be a good indicator of decline 
rate. 



4.5. 4570 Emission Feature 

Similar to Figure [TTJ, Figure [T2] shows what appears to be weak dependence on Am^ 
with the 4570 A emission peak. As with Figures 191 and [TOl the more deviant SNe appear on 
the outer edges of the distribution. This tendency is somewhat stronger than that found for 
the 5150 A feature and is in the opposite direction: the X value for 4570 A is increasing 
with increasing Am 15 . 

There is wide variation in early time evolution to this feature, so the data comprising 
Figure [12] is restricted to spectra taken within 5 days of maximum. Similar to the 5150 A 
feature, the utility of the 4750 A feature in specifying decline rate is uncertain. 

4.6. Ratio Between the Emission Features at 4570 A and 5150 A 

The suggestion of a correlation between Am 15 and the ratio of these two features appears 
in Figure [T3J though it is also weak, particularly for the fast decliners. It appears that the 
4570 A feature has a stronger effect on this ratio than the 5150 A feature. 

Due to the restriction on the 4570 A feature, the same restriction to spectra within 5 
days of maximum is applied to Figure [T3J This ratio may be a useful parameter for the 
slower decliners but not for fast decliners. 
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Fig. 11. — The correlation of the emission peak at 5150 A and Am 
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Fig. 12. — The correlation of the strength of the emission peak at 4570 A and A 
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Fig. 13. — The correlation of the ratio of the strength of the emission peak at 4570 A and 
5150 A and Amis. 
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4.7. Extinction 

To explore the impact of reddening on the spectral indices, each spectra in the supernova 
sample in Table H~T1 was reddened. Four values for E(B-V) were chosen: -0.25, 0.25, 0.5, and 
1.0. An identical procedure to that described in Section [3] was applied to these reddened 
spectra and spectral indexes were recalculated. 

Figures [H] and [15] are similar to Figures [7] and [BJ but now include the effects of various 
values of E(B-V). Figures [T4l and [T5l together demonstrate that reddening has only a minor 
effect on spectral indices and is, in any event, an effect that can be corrected. 

5. Conclusions 

Efforts to correlate the intrinsic brightness of Type la supernovae with spectral line 
ratios or other methods related to feature strength are considerably complicated by line 
blending and noise. We introduce a comparatively new technique to astronomical image 
processing based on an a trous wavelet decomposition and use it to extract spectral strengths 
of Type la supernovae. In a straightforward manner repeated application of the a trous 
generates successively smoother scales. The lowest scale can be identified with noise; a 
smooth residual results from truncating the algorithm at the highest scale. The intermediate 
scales are those that can be identified with spectral features and combining several of these 
intermediate scales provides a robust measure of spectral strength without having to wrestle 
with integration limits or a definition of the continuum. Monte Carlo methods in conjunction 
with a very high SNR Type la supernova spectrum allows for correction of the spectral indices 
of supernovae with lower SNR, even those whose SNR approach one. These same methods 
readily permit error bars to be assigned to the spectral indices. The result is a definition 
of spectral line strength that is applied in this paper to the temporal evolution of spectral 
line strengths and the correlation of important spectral features with Amis. These indices 
are also shown to be largely impervious to reddening. A more robust spectral line strength 
like that developed here is likely to advance the identification or study of Type la subtypes 
or permit the construction of Type la templates with stretches that differ from unity. The 
latter effort, in particular, is directly related to our work on photometric redshift estimation 
for the next generation of ground- and spaced-based survey telescopes. 
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Fig. 15. — Same as Figure El but with the addition of four panels each displaying the 
spectral indices for four values of E(B-V). The lower right panel overlays data from the first 
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